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Diatomic molecule as a quantum entanglement switch 
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.Abstract 

We investigate a pair entanglement of electrons in diatomic molecule, modeled as a correlated double quantum dot attached to 
the leads. The low-temperature properties are derived from the ground state obtained by utilizing the Rejec-Ramsak variational 
technique within the framework of EDABI method, which combines exact diagonalization with ab initio calculations. The results 
show, that single-particle basis renormalization modifies the entanglement-switch effectiveness significantly. We also found the 
entanglement signature of a competition between an extended Kondo and singlet phases. 
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. Quantum entanglement, as one of the most intriguing 
features of quantum mechanics, have spurred a great deal 
of scientific activity during the last decade, mainly because 
it is regarded as a valuable resource in quantum commu- 
nication and information processing [T]. The question on 
entanglement between microscopic degrees of freedom in a 
condensed phase have been raised recently [2], in hope to 
shed new lights on the physics of quantum phase transitions 
and quantum coherence [3] . In the field of quantum elec- 
tronics, a pair entanglement appeared to be a convenient 
tool to characterize the nature of transport through quan- 
tum dot, since its vanish when the system is in a Kondo 
regime [3]. The analogical behavior was observed for two 
qubits in double quantum dot, for either serial and paral- 
lel configuration [S] . The latter case is intriguing, since the 
concurrence |H] at T = changes abruptly from C « f to 
C = when varying the interdot coupling, so a finite An- 
derson system shows a true quantum phase transition. 

Here we consider a nanoscale version of such an entangle- 
ment switch, inspired by conductance measurements for a 
single hydrogen molecule [7] . A special attention is payed to 
electron-correlation effects, in particular the wave- function 
renormalization Recent experiment [9] shows the cur- 
rent through a molecule is carried by a single conductance 
channel, so serial configuration shown in Fig.[T]seems to be 
the realistic one. The Hamiltonian of the system is 
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Fig. 1. Diatomic molecule modeled as a double quantum dot attached 
serially to the leads. A cross-section of the single-particle potential 
along the main system axis is shown schematically. 

H = Hl + Vl+Hc + Vr + Hr, (1) 

where He models the central region, ^^l(_r) describes the 
left (right) lead, and V/^(i?,) is the coupling between the lead 
and the central region. Both i?L(fl) and V^^n) terms have a 
tight-binding form, with the chemical potential in leads /j,, 
the hopping i, and the tunneling amplitude V , as depicted 
schematically in Fig. [T] The central-region Hamiltonian 



Hr 



+ [Zef/R (2) 



(with i,j = 1,2 and a =T,i) describes a double quan- 
tum dot with electron-electron interaction, tij and Uij are 
single-particle and interaction elements, the last term de- 
scribes the Coulomb repulsion of the two ions at the dis- 
tance R. Here we put Z = \ and calculate all the param- 
eters tij, Uij as the Slater integrals [T^ for Is-like hydro- 
genic orbitals 5'is(r) \/ o? ji^ exp(— a|r|), where is 
the orbital size (c/. Fig.[T]). The parameter a is optimized 
to get a minimal ground-state energy for whole the system 
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Fig. 2. Entanglement and transport through the system in Fig.lTlas 
a function of the chemical potential ^l {top panel) and the average 
filling (ni+n2) [bottom panel). Tick (thin) solid and dashed lines 
shows the concurrence C (conductance G) for T = t/9 and t/4, 
respectively. The limits F — > are depicted with dotted lines in the 
bottom panel. The interatomic distance is R = 1.5a. 
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Fig. 3. Concurrence {tick lines) and conductance {thin lines) at 
the half-filled sector (ni+n2) = 2 as a function of the interatomic 
distance R. The remaining parameters are the same as in Fig. [2] 

described by the Hamiltonian ([1]). Thus, following the idea 
of EDABI method [8], we reduce the number of physical 
parameters of the problem to just a three: the interatomic 
distance R, the lead- molecule hybridization F = V'^ /t, and 
the chemical potential fi (we put the lead hopping t = 
1 Ry = 13.6 eV to work in the wide-bandwidth limit). 

The entanglement between electrons placed on two atoms 
can be characterized by the charge concurrence 



C = 2 max i 0, | {cj^Cja) \ ~ \J {nianja){nianja) \ (3) ^ 



where n jo- = 1 — n j^ ■ We also discuss the conductivity cal- 
culated from the formula G = Gq sin^ {E+ - E_ ) /UN [11, 
where Go = 2e^/?i, and E± are the ground-state ener- 
gies of the system with periodic and antiperiodic boundary 
conditions, respectively. Either the energies E± or correla- 



tion functions in Eq. ([3]) are calculated within the Rejec- 
Ramsak variational method 1 1 IJ , complemented by the or- 
bital size optimization, as mentioned above. We use up to 
N = 10* sites to reach the convergence. 

In Fig. [2] we show the concurrence and conductance for 
R — 1.5ao (where ao is the Bohr radius) and two values 
of the hybridization T — t/9 and t/A. The conductance 
spectrum asymmetry, caused by wave-function renormal- 
ization [8] , is followed by an analogical effect on entangle- 
ment, which changes significantly faster for the upper con- 
duction band, where the average filling is (rii-|-n2) « 3 
(one extra electron). The asymmetry vanish when analyz- 
ing the system properties as a function of {ni+n2), show- 
ing it originates from varying charge compressibility Xc = 
d{ni + n2)/dfi « 2/(C/ii -I- U12) ~ l/a. We also note the 
convergence of discussed quantities with F — > to C w 
1 - Kni+na) - 2|/2 and G w Gasm'^{Tr{ni+n2)/2). 

Entanglement evolution with R is illustrated in Fig. [31 
where we focus on the charge neutral section {ni + 712) = 
2. The abrupt entanglement drop follows the sharp con- 
ductance peak for F = i/9, which is associated with the 
competition between double Kondo and spin/charge sin- 
glet phases fT^. For F = t/4 both C and G dependence on R 
become smooth, but the switching behavior is still present. 
Earlier, we have shown that F = t/A is large enough to 
cause molecule instability and therefore may allow the in- 
dividual atom manipulation 8^ . 

In conclusion, we analyzed a pair entanglement of elec- 
trons in diatomic molecule attached serially to the leads. 
Entanglement evolution with the chemical potential speeds 
up remarkably for the negatively charged system, due to 
electron correlation effects. The switching behavior was also 
observed when changing the interatomic distance. 
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